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This writeup highlights the infrared unsafety of the 'midpoint' cone jet-algorithm and provides 
a brief overview of why this is a serious issue. It then shows how one can build a safe (seedless) 
cone algorithm and discusses the potential impact on measurements. 



Two broad classes of jet algorithm are in widespread use at modern colliders: sequential 
recombination type algorithms, such as and Cambridge/ AachenP and cone- type onesPThe 
former take a bottom-up approach to the problem of defining jets, repeatedly combining particles 
that are closest in some distance measure. They work because the proximity measures used are 
closely related with QCD divergences for particle production, and they are much appreciated 
in the e^e~ and ep communities, both for their simplicity and their modest hadronisation 
corrections. Cone type algorithms take a top-down approach, finding coarse regions of energy 
flow (cones) and identifying them as jets. They work because QCD only modifies the energy flow 
on small scales, and so far they have been the preferred type of algorithm in the pp community, 
because of the greater geometrical regularity of the resulting jets and their sometimes lower 
sensitivity to certain components of the non-perturbative underlying event and pileup. 

Cone algorithms have been in use since the early 1980'sl^ and in the early 1990's awareness 
developed^ of the importance for cone-algorithm formulations to satisfy a certain basic set 
of requirements: they must be fully defined, practical in both experimental and theoretical 
contexts, and cross sections must be finite at any order of perturbation theory, i.e. the algorithm 
must be infrared and collinear (IRC) safe. 

Modern cone algorithms involve two main steps: a procedure to find 'stable cones' (a cone 
pointing in the same direction as the momentum of its contents) and a 'split-merge' procedure 
to convert those cones into jets, resolving the problem of cones that have particles in common 
(i.e. that 'overlap'). 
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The most delicate issue with cone algorithms has been that of finding the stable cones. 
A standard procedure had long been to use all particles (possibly above a seed threshold) as 
directions of trial cones, then for each trial cone to use the momentum of its contents as a 
new trial direction, iterating until stable directions are obtained. The drawback of iterative 
procedures is that new stable cones (and jets) may be found if an extra starting point is added. 
This was known to happen with the addition of soft particles in straightforward iterative stable- 
cone searchesP but it had been thought that a trick of adding extra starting points, at the 
midpoints between the cones already found, would lead to a final set of stable cones that was 
insensitive to the addition of extra seeds. Accordingly a recommendation was madJ^I for the 
Tevatron experiments to use such a 'midpoint' iterative cone algorithm. 

It turns out that while the midpoint 
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Figure 1: A configuration (left) in wliicli the midpoint algo- 
rithm {R = 1) gives different jets if a soft particle is added 

(right). 



fix resolves the infrared problems for 
events with two neighbouring hard par- 
ticles, those problems reappear for t hree 
or more neighbouring hard particles'^'^ 
This is illustrated in figH] where in the 
left-hand particle configuration two sta- 
ble cones (and so two jets) are found 
with the midpoint cone algorithm. If 
a soft, ~ IGeV, particle is now added 
(right), it provides an extra seed leading 
to a third (overlapping) stable cone being 
found, and all the cones are then merged 
into a single jet (for / = 0.5). 

The sensitivity to the set of seeds means that the midpoint cone algorithm is either infrared 
unsafe (without a seed threshold) or collinear unsafe (with a seed threshold). This is a serious 
issue, for many reasons: 1) it defeats the purpose of using a jet algorithm in the first place: a 
jet algorithm is supposed to provide a correspondence between the complex hadron level and a 
simple few-parton picture of an event — this correspondence is meaningless if a random 1 GeV 
non-perturbative particle changes the multi-hundred GeV jets. 2) IRC unsafety invalidates the 
theorems that ensure the finiteness of perturbative QCD calculations, because the jets found in 
(divergent, supposedly cancelling) real and virtual diagrams differ. 3) Pragmatically it limits the 
accuracy with which one can meaningfully predict many observables, as summarised in table [U 
and already programs such as NLO JET^ or allow one to go beyond this order when 

using a safe jet algorithm. Therefore the use of a midpoint algorithm squanders the potential 
for accurate predictions that stems from many years of hard theoretical calculations, and forever 
limits the usefulness of data measured with it. 

A solution to the cone algorithm's problems is to carry out an exhaustive ('seedless') search 
for all stable cones. Since additional soft particles do not change the stability of cones, if one 



Observable 


1st miss cones at 


Last meaningful order 


Inclusive jet cross section 


NNLO 


NLO 


W/Z/H + 1 jet cross section 


NNLO 


NLO 


3 jet cross section 


NLO 


LO 


W/Z/H + 2 jet cross section 


NLO 


LO 


jet masses in 3 jets, W/Z/H + 2 jets 


LO 


none 



Table 1: Summary of the order (a* or a^asw) at which stable cones are missed for various observables with a 
midpoint algorithm, and the corresponding last order that can be meaningfully calculated. (Legacy iterative cone 
algorithms, without midpoint seeds, such as JetClu, fail one order earlier). 




Figure 2: (a) Some initial circular enclosure; (b) moving the circle in a random direction until some enclosed or 
external point touches the edge of the circle; (c) pivoting the circle around the edge point until a second point 
touches the edge; (d) all circles defined by pairs of edge points leading to the same circular enclosure. 



has already found all stable cones adding a soft particle cannot lead to extra stable cones being 
found, and so the IRC safety problem is eliminated. A seedless algorithm had been proposecfH^ 
for perturbative calculations, but since it took time ~ N2^ to find jets among particles (10^^ 
years for = 100), it was unthinkable to use it at hadron or detector level. 

Recently it was observe(iI2] that it can be ad- 
vantageous to relate sequential-recombination jet 
algorithms to problems in computational geometry. 
It turns out that this is true also of cone algorithms, 
for which the exhaustive stable cone search reduces 
to a 2D 'all distinct circular enclosures' problem. 
While apparently not having been considered by 
the computational geometry community, this prob- 
lem is easily solved, essentially by considering all 
circles having a pair of particles on their circum- 
ference, cf. fig. [21 With the aid of further standard 
computational techniques one obtaini^l a seedless 
algorithm that takes O (A^^ In A^) time. Not only 
does this provide a practically usable IR safe cone 
algorithm, but it even scales better at large A^ than 
midpoint algorithms (A^^) and is of similar speed 
to them for the values of A^ ^ 500 - 1000 that wih 
be found at low-luminosity LHC. 

Given the cone algorithm's chequered history 
with IRC safety, it is important to establish, as far 
as possible, that there are no further unpleasant 
surprises waiting to be discovered in a few years' 
time. This has been done in two ways: with a 
detailed analytical proof, and via Monte Carlo tests 

in which one finds jets in a 'hard event' (with up to 10 hard particles), repeatedly adds infinitely 
soft particles and verifies that the jets found are the same. If they are not, then the algorithm 
is IR unsafe. The failure rates on this test are shown for a variety of cone algorithms in fig. [3l 
Among the discoveries made in these tests, was that the split-merge procedure also had the 
potential to create IR safety problems. The final version of the seedless algorithm, named 
SISCone, has passed several billion hard event tests without failure. The code for the algorithm 
is available publicl}^^ both in standalone form and as a Fast Jet^l^ plu gm. 

The physics impact of switching from the midpoint to SISCone depends on the observable 
and is illustrated in fig. [H for two cases. For inclusive quantities, like the inclusive jet spectrum 
(upper panel), one sees effects of the order of a couple of percent, as is to be expected since 



1 1 r-| 1 1 r-| 1 1 r-| 1 1 r-| 1 1 — r 


JetClu 


50.1%" 




^^^H SearchCone 


48.2°/d 




MidPoint 


16.4% 




^^^^K Midpoint-3 


15.6% 






9.3% 


Seedless [SM-p,] 1.6% 








^■■■D.17% Seedless [SM-MIP] 


< 10"® Seedless (SISCone) 





10"- 



10"^ 10"-^ 10"' 
Fraction of hard events failing IR safety test 



Figure 3: Failure rates for the IR safety tests with 
various algorithms, including a midpoint variant 
with 3- way midpoints and some seedless algorithms 
with commonly used, but improper, split-merge 
procedures. 
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stable cones are only missed at NNLO onwards. One notes nevertheless that differences of up to 
5% arise when including underlying event effects, and this is related to SISCone's substantially 
lower sensitivity to diffuse 'noise' in an event. 

For more exclusive quantities the differ- 
ences between midpoint and SIS Cone are 
more significant. For jet-mass spectra in 
three-jet events (lower panel of fig-H]), the dif- 
ference starts are LO, and this can translate 
to 40% effects in partonic predictions (which 
essentially corresponds to an unavoidable 40% 
non-perturbative ambiguity for the midpoint 
algorithm) . 

To conclude, while both sequential recom- 
bination and cone- type jet algorithms have 
their place at hadron colliders, it is essential 
that they be practical and safely defined. The 
widespread 'midpoint' cone algorithm is not 
infrared safe, and therefore there are strong 
reasons for discontinuing its use in favour of 
a seedless cone algorithm such as SISCone, 
which is both infrared safe and practical at 
parton, hadron and detector levels. 
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Figure 4: Top: relative difFerence between the midpoint 
and SISCone inclusive jet pt spectra at the Tevatron. 
Bottom: relative diflterence between the midpoint and 
SISCone jet mass spectra, in 3-jet events for which the 
second and third hardest jets are in a common neigh- 
bourhood. 
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